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Abstract 

We derive an exact analytical expression for the one-loop partition function of a string 
in AdS^ X background with world-surface ending on two anti-parallel lines. All quantum 
fluctuations are shown to be governed by integrable, single-gap Lame operators. The first 
strong coupling correction to the quark-antiquark potential, as defined in A/" = 4 SYM, is 
derived as the sum of known mathematical constants and a one-dimensional integral rep- 
resentation. Its full numerical value can be given with arbitrary precision and confirms a 
previous result. 
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1 Overview 

The chance of studying weakly-coupled string theory to gain insight into strongly-coupled gauge 
theory, provided by the AdS/CFT correspondence, has a seminal example in the string realiza- 
tion of the quark-antiquark Wilson loop [1, 2], with heavy quarks modeled by W-bosons. The 
expectation value of the rectangular loop with length T and width L, which in the limit T ^ L 
can be seen as a pair of anti-parallel lines (the "quark" trajectories) at distance L, is given by the 
effective energy of a string on AdS^ x whose ends, restricted to the four-dimensional boundary 
of AdS^, are at a distance L apart. In this context, the potential exhibits a Coulomb-like law 
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where c(A) is a function of the string tension (or 't Hooft coupling) that behaves as 
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Above, the weak-coupling expansion is the field-theoretical calculation of [3, 4], and IK = ^(2) 
is the complete elliptic integral of the first kind with modulus k = A= The problem of finding 

v2 

the first quantum string correction ai to the classical result of [1][2], initiated in [5, 6, 7], has 
been first addressed in [8, 9], resulting in a formal expression for the one- loop contribution to 
the effective action as a ratio of determinants of two-dimensional generalized Laplace operators. 
A numerical prediction for oi has been presented in [10]. Our main motivation here is to address 
the issue of exploiting exact analytical methods for computing the determinants in the partition 
function of [8, 9] and thus the analytically exact value of the constant oi in the sub-leading 
correction to the potential. 

The evaluation of quantum corrections to the energies of classical string solutions in AdS^ x 
[11], crucial device for checking the detailed structure as well as the integrability of the 
AdS/CFT system [12, 13], is in general a hard mathematical problem. The task is simplified 
considering scaling limits of some "semiclassical parameters", as in the case of fluctuations 
over the open string solution dual to the cusp Wilson loop [14, 15, 16], or the closed string 
solutions of [17, 18, 19, 20]. In these limits the solutions become linear in the world-sheet 
coordinates (t, a), thus making constant the coefficients in the fluctuation Lagrangian. In the 
case of the Wilson loop of a pair of anti-parallel lines, which has no other parameters than the 
distance between the lines, the complicated cr-dependence of the lagrangean coefficients makes 
non-trivial the evaluation of the operator spectra. The same is true for the straight line and 
circular Wilson loops [8, 9, 21], for which a first explicit computation of fluctuation determinants 
has been carried out in [22]. There, based on the effective one-dimensionality of the spectral 
problem, it was possible to trade the explicit evaluation of the eigenvalue spectrum for the 
relevant operator with the resolution of the associated differential equation, an approach known 
as Gelfand-Yaglom method [23, 24]. In an analogous fashion the case of the anti-parallel lines 
has been studied in [10], where each functional determinant has been formally expressed in terms 
of the associated initial value problem with Dirichlet boundary conditions (the appropriate ones 
in this framework [25]). While the possibility of a completely analytical treatment of such 
initial value problem was not recognized in [10], the coefficient oi in (1.2) was worked out by the 
authors with great numerical precision. A step forward in the exact analytical treatment of string 
quantum corrections has been made in [26] for the case of the folded string solution [27, 17], and 
recently in [28] for the case of pulsating strings. It has been there realized that fiuctuations on 
this basic class of elliptic solutions can be put into the standard (single-gap) Lame form, which 
allows an exact treatment of the fiuctuation problem. This is useful to extract information in 

^See Appendix B for notation. We adopt here the Abramowitz-Mathematica notation for the modulus of the 
elliptic functions. 



2 



the meaningful semiclassical limits of large [29, 30] and short values of the conserved charges 
typical of the problem [31, 32, 33]. 

We revisit here the evaluation of the one-loop partition function that defines the first sub- 
leading correction to the quark-antiquark AdS potential, showing that also in this case the 
fluctuations are governed by Lame operators. This allows us to present some analytically exact 
results, as the expressions for the fluctuation determinants (2.10)-(2.12) with (2.7), (2.13)-(2.14) 
and the resulting formula for the one-loop partition function (3.2)-(3.3) with (3.4)-(3.8), see 
also the equivalent expressions collected in Appendix D. Finally, we find for ai the following 
representation 

ai = ^-31n2 + ^(]K-^/2(^ + ln2)+X'^"°^) (1.3) 
= -1.33459530528060077364... , 

where the contribution X""™, whose one-dimensional integral representation is displayed in 
(D.IO), can be evaluated with arbitrary precision. The numerical value of ai confirms the 
result obtained in [10]. 

The connection of the fluctuation problem to the integrable Lame differential equation is 
not surprising, since the minimal surface corresponding to the Wilson loop of anti-parallel lines 
belongs, as the folded and pulsating string cases, to the important class of classical string 
solutions expressed in terms of elliptic functions (see Appendix A). It is however interesting to 
see on this non-trivial example how the integrability of the a-model on AdS^ x [34] [35, 36] is 
extended from the classical to the one-loop level via this special, integrable, type of potential. 
It is also interesting to recall that the chance of exploiting the integrability of the underlying 
sigma-model to calculate Wilson loops within the AdS/CFT correspondence [37] has been made 
recently concrete, via the connection of Wilson loops to A/" = 4 SYM scattering amplitudes, with 
the proposal of Thermodynamic Bethe Ansatz equations for the latter [38]. Although the Wilson 
loops appearing in amplitude computations consist of light-like segments that are not obviously 
related to the configuration of space-like anti-parallel lines of interest here, it is fascinating to 
think about the possibility (on the lines of the approximation in [39]) of using a description 
similar to the one in [38] also in this case. 

It would be also interesting to exploit similar analytical methods in the case of the one-loop 
partition function for the anti-parallel lines configuration in a Schwarzschild-74(iS'5 background [6, 
40], whose formal expression has been worked out in [41]. 

The main body of this article contains the analytical study of the fiuctuations and the evalua- 
tion of the one-loop contribution to the quark-antiquark potential. Appendices A, B and C recall 
basic facts on the world-sheet set-up, elliptic integrals and functions and on the Gelfand-Yaglom 
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method. Alternative expressions for the relevant integrals are displayed in Appendix D. 



2 Fluctuation operators and their analytical determinants 

Given the invariance of the anti-parallel lines configuration under time-translation, both the 
bosonic and the fermionic fluctuation lagrangeans depend non-trivially only on the cr-coordinate, 
and the original two-dimensional spectral problem is reduced to the evaluation of one- dimensional 
functional determinants. After suitable world-sheet reparametrization and fermion diagonaliza- 
tion [9, 10] reviewed in Appendix A, the resulting effective action for a string in AdS^ x 
background with world-surface ending on two anti-parallel lines [8, 9] can be written as follows 
upon Fourier transformation of the time variable {dr = —i oj) 

f duj det^O+det^O. 

where T = / dr is the r-period. Above, Oq = —d"^ + uP' is the free operator and 

O, = -dl\V,{a) + 0? i = l,2,±, (2.2) 



1 ^^1 2 T. sno" dncj 

, Vi = Cn £7, V± = — 2 

cn^fj cn^cj Zenker 



The Jacobi elliptic functions appearing in (2.3) and defined in Appendix B have fixed modulus 
/c = and — K. < cr < K. The operators Cj are also defined in (A.10)-(A.14). 

The partition function (2.1) suffers in general from linear infrared divergencies, that can be 
cured subtracting a reference solution (as in [22] for the example of the circular Wilson loop). 
The one-loop correction to the quark-antiquark potential can be thus obtained subtracting twice 
the infinite, self-energy contribution of each of the parallel lines (quarks) [10], and dividing over 
the infinite time period T = J dt 



V„^P = lim ^ 



2r| 



■^■r^oo, (2.4) 



TT 

3 



where the relation between T and T follows from (A. 5) and (A. 8) 

Exploiting elementary transformations of Jacobian elliptic functions [42], it is easy to check 
that each non-trivial fluctuation operator is a single-gap Lame operator with the following eigen- 
value problem 



dt + sn^x + 



/a(x) = A/a(x) , (2.5) 



''The world-sheet of the straight hne can be parametrized with the same time variable as the one for the 
anti-parallel lines [10], see (A. 8). 
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where, with respect to (2.2), x and are a shifted (and rescaled) a variable and euchdean 
frequency respectively. Explicitly, 



(a) for the first bosonic operator Vi: x = a + {1 + i)'K and il^ = — 1; 

(b) for the second bosonic operator V2: 2; = (1 + i) o" + K and = ^ — 1; 



(c) for the fermionic operators : x 



, and S Z = 2 1 w — 1 

[f(-l + i) + f(3-i), fory+ 



The Lame spectral problem (2.5) can be solved exactly, and hence the corresponding determinant 
can be computed analytically, relying on the knowledge of the solutions to (2.5) and the use of 
the Gelfand-Yaglom method (see, for example, [24] and [26] [28]). While the general procedure 
is briefly reviewed in Appendix C, let us see explicitly the evaluation of the determinant for the 
bosonic fluctuation governed by the first potential Vi. 

Readapting the solutions (C.5) to the case (a) above, two independent solutions of the relevant 
differential equation are 

(„\ = H{a + K{l + i)±ai) ^:^z{ai) {a+K(l+t)) ^ ^3(^^''2k''^) gT^(ai) (<T+K(l+i)) (2 6) 

e{a + K{l + i)) ^2(11) 

where the Jacobi H, and Z functions are defined in (B.8) in terms of the Jacobi ^-functions 
and 



«i = sn-^ = sn' Vl + 2lo^ . (2.7) 

The solutions (2.6) diverge at the extrema a = — IfC and a = +IK of the interval, which is a direct 
way to see the standard need [25] of an infrared regulator e. The Gelfand-Yaglom theorem will 
be therefore applied solving the initial value problem in the interval — IK + e<(T<IK — e where 
e is arbitrary small. The linear combination 

f ^ _ y+i-^ + e) y-{x) - y-{-^ + e) y+jx) 

^ W{-^ + e) ' ^ 

with wronskian W{x), see (C.7), evaluated at the regularized initial point, is a solution of the 
homogeneous equation with boundary conditions 

ti(-]K + e) = 0, u'{-'K + e) = l . (2.9) 

As follows from the discussion in Appendix C, the determinant of the bosonic operator Vi with 
Dirichlet boundary conditions in the interval [— K + e, K + e] will be then given by u{K — e). One 
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finds 



^ ^'"^'^ "^'"1 sinh[ 2 (K - e) + Si ], Si = In . (2.10) 

snaicnoidnai L V iM J ^ ^ 



In a similar fasliion one can work out the regularized determinants for the potential V2 and 
V±, obtaining 

1 2ns2[(l + i)e] -sn2a2 



detOo 



sinh 



ITT a2 



(1 + i) sna2 cna2 dna2 

sna/cna/dna/ ^iImTTi) 



2Z(a2)(l + i)(]K-e) + — ^ + S2 , (2.11 



where 



(IfC — e) ITT af 



+ i 



+ S/ ,(2.12) 



02 = sn ^\/l — 



a J = sn ^ \/l + 4iu/^ , 



and 



In- 



^ 7r(a2 + (l+»)£) \ 
k 2K / 

^ 7r(a2-(l+t)e) -j ' 



2K 



In 



^4(^(«/ + T^)) 

^2(m(«/-tti))) ■ 



(2.13) 



(2.14) 



Notice that in the fermionic case det 0+ = det O- = det Of. This can be understood by noticing, 
in (2.3), that V+{-a) = VI (cr). Namely, it holds that 0+ = V'^ O- V, with V the unitary parity 
operator with respect to a, implying det 0+= det O- ^. 

The contribution of the massless bosons can be easily obtained via the same method 

sinh[2a;(]K - e)] 



detOn 



(2.15) 



Expanding in e ~ and retaining the divergent contributions, one gets 

2 



detC^ 



detO) ^ 



detO^ ^ 



e2t^2 V 4a;4 - 1 
1 



sinh[2KZ(ai)] , 



£2^2 V W^ + l 

4 1 

sm 



sinh 



2(l + i)KZ(a2) + 



^ TT 02 



smh[2Kuj] 



(l + i)]KZ(a/) + 



vra/ 



(2.16) 

(2.17) 
(2.18) 

(2.19) 



''in this case, in which the evaluation of determinants is done via the Gelfand-Yaglom theorem, the equivalence 
of the determinants can be easily checked exploiting this parity feature as inherited in the solutions via which the 
determinant is defined, see (2.8) and (C.7). 
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As usual, the divergence ~ 1/e in the resulting ratio of determinants as they appear in (2.1) 
is cured subtracting twice the contribution of the straight line, which can be evaluated and 
regularized by the same means [22] [10] 



r 



r 

' 2 



2^ 



In 



1 + 



1 



UJ € 



(2.20) 



3 One-loop correction to the quark-ant iquark potential 

The one-loop correction to the quark-antiquark potential is formally defined by (2.4), in which Fn 
is given in terms of the determinants (2.16)-(2.19), is substituted by the regularized expression 
(2.20) and the regulator e is sent to zero. Namely, it is 



(1) 



lim 

2 KL e->0 









In 


Loo 2^ 





det^O 



/ 



2 In 



1 

1 + — 

UJ € 



(3.1) 



(detOf)2detO| det^Og- 

Making (3.1) explicit in terms of the determinants (2.16)-(2.19), the following exact analytical 
expression is obtained for the one-loop correction to the quark-antiquark AdS potential 

^ r+oo 

2KlJo 



(1) 



qq 



duj In 



1282^10(1-4^4)^/1+^ 



(16^4 + 1) 



+ 



(3.2) 



1 



2KL 



vr 



duj In 



■I 



sin^ 


\l + i)KZ{af) +^ 




sinh2[2]fCZ(ai)] sinh 


2(l + i)]KZ(a2) + ^" 


sinh^[2Ka;] 



,/2KL 2KL ■ ^^'^^ 
Above, ai,a2,af are implicitly defined in (2.7), (2.13) and in the last equivalence we have 
reported the result for the integral in the first line. The second non trivial integral I can be 
partially given in terms of known mathematical constants. One can proceed rewriting it as ^ 



X =1 doj In 
/o 



cosh® X f 



■oo 

+ / duj In 
'k 



cosh X 



f 



sinh xiinf sinhx2 sinh [2]Kl(j] 



-cos^xizcro sinhx2 sinh^[2]Ka;] 

where 

xi,ero = K (f[cos-1(V2u;)] - 2E[cos-\V2u;)] ) + ^F[cos-H^g)] , 
xiir,i = K(F[sec-\V2ij)] - 2E[sec'\V2 
F[cos"i(i-^' 



(3.4) 



X2 



:uj) 

,2 , 



2i?[cos-Hi-S)]+^^^ 



+ lkn^os~H'T^)] , (3.5) 
)+f-^F[sec-n^)] , (3.7) 



''A quick way to obtain tiiis expression consists in taking the derivative of the arguments of the hyperbolic 
functions in (3.3) and integrating back using standard tables of integrals. 
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and E[x] = E[x, ^] and F[x] = F[x, ^] are the incomplete elliptic integrals of the first kind 
defined in (B.l). The need of two separate intervals of integration, u £ {0,k) and uj G {k,oo), 
is due to the bosonic fluctuations described by Oi and is clear, for example, by looking at the 
square root in (2.16). 
Rewriting now 

Insinhx = In | + x + ln[l - 6"^^], (3.9) 

its analogue for the cosh, and disregarding the constant contribution which will vanish due to the 
balance of the world-sheet degrees of freedom, one can consider the part ~ x in (3.9) and work 
out some analytical finite contribution (the logarithmical and power-like divergencies will cancel 
in the ratio). The numerical integration for the remaining contribution ~ ln[l — e~^^] converges 
quickly to a steady value, and can be obtained via standard packages like Mathematica with an 
arbitrary precision. This way the contribution of each fluctuation can be evaluated separately. 
For example, in the case of the fermions the indefinite integration will give 

8 /" duj'xf = 8uj'^K- — lnuj + 2K--(2 + 3ln2) + o(-^). (3.10) 

Jo JK JK \uj^J 

In an analogous way one can evaluate the analytical contributions for all the fluctuations, check 
the cancellation of the divergent pieces and get for the finite ones 

Werm = 2IC - ^ (2 + 3 In 2) , Xf^^^ = ^, (3.11) 

Xr = J(l + ln2)-V2 1n2, = ^ (i + i^) _ K , (3.12) 

where Ifree amounts for the total contribution of the free fiuctuations to I. 

The remaining contributions can be evaluated numerically with arbitrary precision ^ 

1.41586 , (3.13) 

/•oo 

2 dhj ln[l - exp(-2xiinf)] = 1.18174 , (3.14) 

Jk 

-- 0.43859 . (3.15) 
Adding together the analytical and numerical contributions 

X^'^ = -f— --ln2) +lC-^ln2, X'^"'" = 3.09111 , (3.16) 



T-num 

-^ferm 


poo 

8 duj ln[l + e-2^/] 


Jo 

rk 


T-num 


— 2 I duj ln[cOS Xizoro] 
JO 
/■oo 


-7-num 

-^2 — 


- dw ln[l 
Jo 



K V 24 2 



^Notice that the first term in If"™ automatically includes the type of constant contribution, In i , which should 
cancel in the balance of degrees of freedom. Such contribution, amounting to 2 dojln | = — In 2, has then 
to be subtracted and is in fact included in If". 
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it follows for the one-loop correction in (3.3) the expression 
.(1) _ 1 



'vr /57r 








.kI 24 





^/2{^T + In 2) + X" 



0.30492 
Z 



(3.17) 



2KL 

A compact way to define X"^*^™, which is equivalent to the sum of the contributions in (3.13)- 
(3.15), is given in (D.IO) and evaluated there with high precision. The AdS quark-antiquark 
potential is therefore given by 



1 



A > 1, 



(3.18) 



where the one-loop correction oi is given in (1.3), and confirms the result obtained in [10]. As 
observed there, it is interesting to notice that, when compared to the strong coupling prediction 
via summation of ladder diagrams of [3] 

1 



vr L 



qq 



Va 



A > 1, 



(3.19) 



ai has the same sign and smaller absolute value. 

To summarize, rephrasing the fluctuations over the minimal surface related to Wilson loop 
with anti-parallel lines in terms of the Lame spectral problem (2.5), we were able to present useful 
analytical formulas for the fluctuation determinants (2.10)-(2.12), for the partition function 
(3.2)-(3.3), and finally the representation (1.3) for the first subleading correction to the quark- 
antiquark potential. 
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Appendix A: World-sheet set-up 

At the classical level, the quark-antiquark potential in AdS^ x described by the metric 

ds^ = y^{dx^f + '^ + dVL\ , (A.l) 

9 



is evaluated considering two anti-parallel lines extended in the direction and located at = 
ib-^. In (A.l), we set to 1 the radius of both AdS^ and S^, n = 0, 1, 2, 3 and the index 4 labels 
the coordinate y, ranging from its minimal value in the bulk to an infinite value on the boundary. 
The world-sheet in the bulk can be parametrized by (r, o") = {x^,x^), —T/2 < r < T/2. Given 
the invariance of the problem under x'^-translation, the surface stretched between the lines is a 

function y = y{x^) of the x^ coordinate only, and the induced metric reads 

1 

? 

Given the Nambu-Goto action 

Vat 



ds^ = y^dt^ + ^{y^ + y'^)da\ (A.2) 



S 



I daVy'^ + , (A.3) 



2tt 

its equation of motion y y" = Ay'"^ + 2y^ has a first integral of motion 

y'^ = 4 - y\ (A.4) 

which can be can be integrated in terms of elliptic functions. Above, y^ is an integration constant 
corresponding to the minimal value of the coordinate y in the bulk and is related to the distance 
L between the lines via 

One then proceeds evaluating the action (A.3) on the solution (A.4) 

day — > S = -——y-, (A. 6) 



2vry2 7_^/2 41K2 L' 

where, following [9], one notices that (y'^y')' is a total derivative and replaces y'^ by —yQ, thus 
assuming that the infinite boundary contribution can be dropped. Such a prescription coincides 
with normalizing the partition function to the straight line case [25]. The classical contribution 
to the quark-antiquark potential is obtained dividing the action by the infinite time period 
T — )• oo, thus obtaining Vqq{L) = — j^t; as in the leading part of (3.18). 

The one-loop correction to the result (A. 6) is obtained by considering fluctuations over the 
classical solution, a problem addressed in [8, 9]. Bosonic fluctuations are obtained via a standard 
background field method, while a cj-dependent rotation in the target space and the standard k- 
symmetry gauge fixing 9^ = (P' for the two Green-Schwarz spinors are used to put the quadratic 
fermionic term in the Green-Schwarz action into the standard kinetic term for a set of 2-d 
Major ana fermions. In static gauge ^ the resulting one- loop partition function is 

detS/2(_^^ay^ + ^g) 



det2/2(_v2 + 2) det^/2(-V2 + \m + 4) det5/2(-V2) 



(A.7) 



^The infinite contribution of the ghost determinant is regularized changing the normalization of the non-trivial 
(longitudinal) bosonic fluctuation, as seen in details in [9]. 
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where R^'^^ is the scalar curvature, 7'^ = T2, 7^ = ti and 7071 = —iT^ are the Pauh matri- 
ces. As suggested in [9], it is useful to deal with a conformally flat induced metric, obtained 
reparametrizing the world-sheet via Jacobi elliptic functions of fixed modulus k = ^ [10] 

y = — > i = . (A-8) 



cncr' -y/2 



2/0 



where now — K < ex < IfC and — ^ < r < ^. 



The induced metric and the scalar curvature read then 

dsfnd = J^(dT^ + ^^^)' ^^^^ = -2(1 + cnV) . (A.9) 

The explicit expressions for the bosonic differential operators appearing in (A. 7) are then [10] 

- = -2cnV {d^ + dl) = 2 cnV Oq (A. 10) 

-V2 + 2 = -2cnV(52 + 9^) + 2 = 2cnVC'i (A.ll) 
-V2 + 4 + i?(2) = -2cnV(c)^ + a^) + 2(l-cnV) = 2cn2(TC'2 (A.12) 

where the operators Oq, Oi and O2 are defined in (2.2)-(2.3) upon Fourier transform of the time 
variable {dr = —iuj). As suggested in [10], the fermionic differential operator 

-«7 Va -hT3 = V2cncT -i{da + — ]ti-ojt2 + ^ rs =^2cTiaO^ (A.13 

L V 2cno- / V2cn(j -I 

can be further diagonalized after squaring it. Using for example M = J one has 

Ol = yHiT^Afdiag{0+, 0-}M^^= , (A. 14) 

^ Vcncr 

where 0+ and 0+ are defined in (2.2) and (2.3). Therefore, it is det^/^O^ = det^C'+ det^O.. 

Each "flat-space" operator O above is rescaled with respect to the original differential operator 
appearing in (A. 7) via the measure = 2cn^a. The finite contribution of such measure to 
the logarithm of the original determinant (related to the Seeley coefficient which determines the 
conformal anomaly ^) can be explicitly shown to cancel in the ratio (A. 7) of determinants [10]. 
This justifies the final expression (2.1) of the effective action. 

Appendix B: Relevant elliptic functions and identities 

The incomplete elliptic integrals of the first and second kind are defined via 

F[x, k^] = I de {I- sin^ 9)-^''^, E[x, k"^] = dO {I - k^ W Of'^ (B.l) 
Jo Jo 

*See discussion in Appendix A of [9]. 
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where fc^ is their modulus. The corresponding complete elhptic integrals are given by 



K{k'^) = K = F[§, , E{k^) = E = E[^,P] . 

Defining the Jacobi amplitude as 

(f = am(ti I k^), where u= dO {I - k'^ sin^ ey^/"^ 

Jo 

the Jacobi elliptic functions sn, cn, dn are defined by 

sn(n I k"^) = sin tp, cn{u | /c^) = cos if, dn(n | k"^) = [1 — k^ sin^ y^Y^'^ 

and, for example, ns(ti | /c^) = l/sn(ti | k'^). 

Useful relations between the squares of the functions are 

-dn2(n I k^) + k'^ = cn\u \ k^) = ^ sn^{u \ k^) - k^ 
-k'^ nd(n I k^) + k'^ = -k^ k'^ sd^iu \ k^) = k^ cd{u \ k^) - k^ . 



(B.2) 



(B.3) 



(B.4) 



(B.5) 
(B.6) 



(B.7) 



A useful identity is 

sn~-^(z, ^) = F(sin~^ z, ^ 
The Jacobi H , and Z functions are defined as follows in terms of the Jacobi 9 functions 

«(„ I,. , = e(„|.^) = ..(||., 

where q = q{k'^) = exp(— vr^^). A useful identity is 



E 

Z{x\k'^) = E{x\k^) F{x\k'^) . 



(B.f 



Appendix C: Lame problem and determinant via Gelfand-Yaglom 
method 

Following [24], consider a Schroedinger operator on the interval x G [0, L] with Dirichlet bound- 
ary conditions 



dl + V{x) ij{x) = \^l){x), ^(0) = 0, ^'(0) = 0. 



(C.l) 



Then to compute the determinant one has to solve the associated homogeneous initial value 
problem 



dl + V{x) ct>{x) = \ct>{x), m = ^. '^'(0) = 1 



(C.2) 
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and 



det 



dl + V{x) =<p{L) . 



(C.3) 



For the single-gap Lame problem 



5^+ 2k'^sn^{x\k^) /(x) = A/(x) 



two independent solutions are [43] 

/^(^■)-^'^'^<°'' V 

In terms of them, a solution satisfying the conditions in (C.2) is 

1 



sn(a I k 



1 + A;2 - A 



u{x; A) 



W{x) 



U{x) Ux) - Ux) U{x) 



where W is the wronskian at a generic initial point x 

W{x) = U{x)f'_{x)-f'+{x)f.{x). 



(C.4) 



(C.5) 



(C.6) 



(C.7) 



Exploiting f±{—x) = — /=p(— x) and some properties of the Jacobi elliptic functions it is then 
easy to check that, in the interval [— IK, K], the expression for the determinant (C.3) yields ^ 



DetDir = u{K; A) 



cna 



sno dna 



sinh[2]KZ(Q)] 



sinh[2ICZ(a)] . (Ci 



- AVl + A:2 - A 

The determinants (2.10)-(2.12) evaluated in Section 2 are generalizations of the expression (C.8). 
Their slightly more involved form is simply due to the presence of the infrared regulator e, which 
alters the boundary conditions of the problem. 



Appendix D: Equivalent form of the integral T in (3.3) 



Basic manipulation of the special functions in (3.3) with identities such as (B.7) and (B.9) leads 
to the following expression 

cosh^ X f 



k 

doj In 







cosh X f 
sin^ xizero sinh5;2 sinh^ [2 IC w] 



•CO 

+ I doj In 



sinh Xiinf sinh5;2 sinh [2]Ka;] 



D.l 



*In the square roots at the second equivalence the known eigenvalues of the Lame equation appear, see for 
example [26]. 
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where 



X2 



H ^ 2 



-u) TP rl 1 5. A a, A 



2' 4' 4' 

3 



-16 a;'' 



^ ^-40;^ 

L 2 ' 4 ' 4 ' 



+ 



8 2i/''l 



5(1+^1-4^4)5/4 



2i^l 



3 5 9.1 

4' 4' 4' 2 

2E 



(1 - Vl - 4^4) 



™ — K/1 _|_,/77T4 2E p ri 1 5. 1 1 I K p f 1 3 7 . 1 1 , 



+ 



K(2a;^+V4a;4„ 1)3/4 



5 V2 (8i^4--l+4V4w4_la;2) 2-^1 
2-t'l I o, 7, 7; — W| + 



3 5 9.1 
4' 4' 4' 2 



V4a;4-1 



4aj^ 



TTCJ 77 ri 1 5. , ,41 I lijfK 77 ri 3 7. ,,41 
2-^1 L2'4'4'~^ J 3~2^1 L2'4'4'" 



(D.2) 

(D.3) 

(D.4) 
(D.5) 



A compact way to write the numerical contribution X"^"™ in (3.16) is obtained as follows. 
Each of the (3.5)-(3.8) can be put in a simpler form with the change of variables uj = \ tan ^ in 
(3.5), uj = ^ cos ^ in (3.6), cj = sec ^ in (3.7) and lo = tan ^ in (3.8). One obtains 



^Izero 
^linf 

^2 



vr 

2K 
vr 

m 

vr 



F[f] + 2IfCi^[f] 



+ K ) F[f ] - 2IK E[f ] + lEC tan f V3 + cos 



a 



+ IK) F[al - 2IfC£;[al +]fC sina Wl + tan4- 
21C / ^ ^ ^ ^ V 2 



(D.6) 
(D.7) 
(D.8) 

(D.9) 



which makes explicit = 2x/. In terms of the variables above, the contributions (3.13)-(3.15) 
are summed as 



da 



2 cos 



2 a 



In. 



fl + p~'^^S\'^ / 

- ^ sinf ln(l - e-2^""f) + cos^ f Insinxiz 



1 - 



3.09111054729005989778296487945453992761532660548813 . (D.IO) 
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